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Abstract 

We study solitons arising in a system describing the interaction of a two-dimensional discrete 
hexagonal lattice with an additional electron field (or, in general, an exciton field). We assume that 
this interaction is electron-phonon-like. In our previous paper [J], we have studied the existence 
of two-dimensional solitons and have found that these solitons exist only if the electron-phonon 
coupling constant is sufficiently large. In this paper, we report the results of our investigation for 
small values of this constant, close to its critical value for the existence of solitons. We find that as 
the coupling decreases the soliton gets very broad and then becomes effectively one-dimensional. 
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I. INTRODUCTION 



Recently, a Frohlich Hamitonian was studied on a two-dimensional, discrete, quadratic 
, respectively hexagonal lattice |4|. It is well known that the interaction of an 
excitation field such as an amide I- vibration in biopolymers or an electron (in the case of 
the Frohlich Hamiltonian) with a lattice whose distortion can be caused by this excitation 
results in the creation of a localised state which, in what follows, we refer to as a soliton. 
Such a soliton was first introduced by Davydov p in the 1970s to explain the dispersion free 
energy transport in biopolymers (see also Q for further details). In 

Baa 

, the existence 

of such solitons was studied numerically and it was found that their properties depended 
crucially on the magnitude of the electro n-phonon coupling constant. In particular, localised 
2-dimensional structures exist only above a critical value of this coupling. These numerical 
results have been confirmed by a simple analytic argument: using a Gaussian in 2 dimensions, 
we can approximate very well this critical value. 

In this paper, we extend the results for the soliton on a hexagonal lattice .4] and look at 
the properties of the solitons close to the critical coupling. We find that close to this critical 
coupling, three types of solutions exist: the "narrow" 2-dimensional solitons, discussed in 
^], "broad" 2-dimensional solitons, which extend over much bigger parts of the lattice than 
the former, and one-dimensional solitons. Since our model is a toy-model for solitons on 
nanotubes, these latter configuartions can be interpreted as ring-like solitonic structures 
around the tube. 

In section II we present our equations and in section III we discuss the results of our 
numerical investigations. 



II. BASIC EQUATIONS 



The equations describing an electron field interacting with a deformable hexagonal lattice 
were derived in |4|. These equations read: 
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for the electron field and 
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for the displacement fields U and V in the x- and, respectively, ^-direction. % and j index 
the points in x-, respectively ^-directions. 

? or the corresponding Hamiltonian and the appropriate rescalings, we refer the reader 
to 4J]. Note that K x is the self-coupling constant of the displacement fields and g is the 



electron-phonon coupling. Wq is the phonon energy, while Eq was chosen to be 0.142312 
in our numerical investigations. We have normalised the electron-function ip to one, i.e. 

N X ,Ny 

J2 ^Piji^lj — 1 thus assuming the excitation to be one excess electron. Note that N x and 
N y denote the number of lattice points in x-, respectively ^-directions. 



III. NUMERICAL RESULTS 

In this paper, we present new results on a model describing the interaction of one excess 

n 

electron with a 2-dimensional hexagonal lattice. This model was studied first in |4|. Before 
we discuss our new results, let us summarize the results of |4j. In j4| it was found that 
2-dimensional solitonic structures which are created due to the interaction between the 
distortions of the hexagonal lattice ( "phonons" ) and the electron field on this lattice ( "excess 
electron"), exist only when the electron-phonon coupling is large enough. This numerical 
result was also found in the case of the quadratic lattice and thus seems to be a 

generic feature of 2-dimensional lattices. For values of the coupling constant smaller than 
some critical value g cr , the interaction between the lattice and the electron (exciton) field 
becomes too small for 2-dimensional solitonic structures to exist. In this paper, we present 
new results on the behaviour of the system close to this critical value g cr . We have confirmed 
that "narrow", i.e. well pronounced 2-dimensional solitons exist for g > 2.296. A typical 
solution of this type is shown in Fig. Q]for N x = 160 and N y = 20. 
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FIG. 1: A typical "narrow" 2-dimensional solitonic structure on a N x x N y = 160 x 20 lattice for 
g close to the critical value g cr = 2.296. 

Clearly the probability density ipip* ls localised only on a small portion of the lattice. For 
g e [2.29475 : 2.32592], a broad two-dimensional soliton appears (see FigJSJ), while decreasing 
the value of g further leads to a configuration in which the soliton becomes effectively one- 
dimensional. A typical one-dimensional configuration is shown in Figj2j 

1 -dimensional soliton 




FIG. 2: A typical one-dimensional solitonic structure on a N x x N y = 160 x 20 lattice for g = 2.295. 
As is obvious from this figure, for small g, the function ip^j becomes effectively inde- 
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pendent of the ^-direction and we see N y copies of the function While for g — > oo, 
we expect the solitonic structure to be effectively fixed to one point on the lattice, thus 
(ipi t j?Plj)max ='■ (V , V'*)mox — ¥ 1> the threshold of our one- dimensional soliton is characterised 
by the value of (ipip*) ma x, which can be obtained for the last possible "broad" soliton di- 
vided by the number N y . In our case N y = 20 and the last value of the broad soliton 
(ipip*)max ~ 2.84 • 10 -2 , so we would expect solitons with heights on the order of 1.42 • 10~ 3 . 
Indeed, for g < 2.29475, we find solitons with heights of ~ 1.62 ■ 10 -3 , which confirms our 
interpretation. We plot the g dependence of the value of (ipip*) ma x m FigED Clearly, we 
note the familiar structure of the curve, especially, we notice that at g ~ 1.0, the value of 
(ijjil>*)max becomes constant and is equal to (ipip*) ma x = 3.125 • 10~ 4 . This is not surprising 
since this is just the value where ipip* = l/(N x N y ) = 1/(160-20) w 3.125- 10" 4 , i.e. the field 
is completely delocalised, i.e. it is equally distributed over the whole grid. 

The fact that the one-soliton ceases to exist is, however, related to the number N y . In 
principle, one-d solitons should exist for all values of 2.29475 > g > 0. Of course, this value 
was derived for ipi jipij = 1- To study the dependence on N y , we have looked at the 
dependence of the critical value of g, g cr ,id(^y) ^ which the one-d soliton becomes equally 
distributed over the grid, i.e. ipijipij = ^-/{N x N y ). The values are listed in Table I. 

TABLE I: Dependence of g cr ,id ° n N y for N x = 160 



N y Ny 1 9cr,ld(Ny) 

20 0.050 1.00 

16 0.062 0.88 

12 0.083 0.67 

6 0.167 0.30 



We notice that the value g cr ,id(N y ) decreases with increasing Ny 1 . Since iV" 1 represents 
the normalisation of the electron field function in the x-direction, i.e. of the one-d soliton 
for one fixed y, this means that the "bigger" the one-d soliton is the lower is the critical 
value of the electron-phonon coupling. We can understand this behaviour by observing 
that by scaling the ip field in our equations CQ)-© such that ipi,jip*j — > N^^ipi^ip*^) we 
obtain the equations for a system with a normalization of the ip fields corresponding to 
the one-d solitons. We notice that this rescaling is nothing else but the rescaling of the 
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FIG. 3: The g dependence of the height {ip r tp*)max of the two-dimensional and one-dimensional 
solitons. Note that there exist "narrow" (dashed line) and "broad" (solid line) two-dimensional 
solitons. The figure on the right presents the zoom of the one-d soliton curve (dotted-dashed). 

coupling constant g — > N^g. Following this argument, the product N^g^xd should thus 
be constant. Looking at Table I, we indeed find that N~ l g 



or, Id 



0.05. 



We have also plotted the corresponding energies of the solutions. These are shown in 
Fig0] Note that the configurations with equally distributed ip would have an energy « 
—5.8576. The energies of the soliton solutions are clearly below this value. The energy 
of a narrow soliton is the lowest energy configuration as long as the coupling constant is 
above the value at which one-d solitons can exist. Then, these one-d solitons have the lowest 
energy. Note that the broad soliton never becomes the lowest energy configuration and so 
appears to be unstable. 

IV. CONCLUSIONS 



In this short note we have reported results of our studies of an electron phonon system on a 
hexagonal lattice which was longer in one direction than the other (i.e. with a larger number 
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FIG. 4: The energy of the two-dimensional narrow, respectively, broad solitons in comparsion with 
the energy of the one-d soliton. Note that the energy of the configuration with equally distributed 
ipi,jip*j = 1/(160 • 20) would have energy —5.8576. 

of lattice points in the x-direction - N x than in the ^/-direction i.e. N y ). We have looked at 
the values of the electron-phonon coupling constant at which the system possesses soliton- 
likc solutions. We have found that ''genuine" two-dimensional solitons. like those studied 
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exist for the coupling constant larger than some critical value. For the values of the 
constant close to this critical value the system possess another solution - which corresponds 
to a "broad" soliton. This new solution appears to be unstable. When the coupling constant 
is below its critical value the two-dimensional soliton is so "broad" that it can be interpreted 
as an effective one-dimensional structure. As the value of the coupling decreases further the 
"one-dimensional" soliton becomes broader and, below some new critical value, it spreads 
over the whole lattice leading to complete delocalisation. 

The value at which this happens depends on the size of the lattice. However, this is not 
surprising as for an effective one-dimensional soliton the normalisation of the electron field 
plays an important role (the soliton is normalised on the two-dimensional lattice and so the 
field of the effective one-dimensional soliton is normalised to where N y describes the 
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number of points in the "shorter" direction). A simple scaling argument then explains the 
observed N y dependence of the critical value of the electron-phonon coupling constant for 
the existence of one-dimensional solitons. 

Though our investigations were performed only for the 2-dimensional hexagonal lattice, 
we would expect that similar features appear in the case of other lattice structures, i.e. 
lattices in physical three dimensional spaces, the quadratic lattices of DDD etc. 
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